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Abstract 

We review recent developments in the application of perturbative QCD to phenomena at small x. 
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Both HI Q and ZEUS g] presented measurements 

(N of F 2 {x 1 Q 2 ) obtained from the 1993 HERA run. A 

^ sample of these data is shown in Fig. 1. The dramatic 

! rise of F 2 (x,Q 2 ) with decreasing x, discovered in the 

£C) 1992 data, is now firmly established. The data are 

ON compatible with Altarelli-Parisi (or GLAP) evolution 

from "starting" parton distributions in which the sea 

\j quarks have the small x behaviour 
ON * 



a: 



xS 



with A ~ 0.3, 



(1) 



i see, for example, the MRS (A) curve H in Fig. 1. 



Also OPAL presented a measurement of the 
c~| photon structure function i 7 ^ at x — few xl0~ 2 . No 
^ rise with decreasing x is seen at this value of x. However 
both HERA and LEP2 have the potential to measure the 
photon structure function at much smaller x. Here we 
^ concentrate on the structure of the proton since many 
results have already been obtained in the HERA small 
x regime (x < 10 -3 ). 
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1. Deep-inelastic map 

To orientate ourselves we show in Fig. 2 a map of 
the kinematic regime for deep inelastic electron-proton 
scattering. 



Altarelli-Parisi (or GLAP) equations which are typically 
of the form 



dg/d\ogQ 2 = P gg ® g + 



(2) 



where the convolution is over the longitudinal momen- 
tum fraction. For simplicity, we concentrate on the 
gluon, the dominant parton at small x. Effectively the 
GLAP equations resum the leading (a s logQ 2 )™ contri- 
butions, which correspond (in a physical gauge) to the n- 
rung gluon ladder of Fig. 3. In fact the leading log arises 
from the strongly-ordered region of transverse momenta 



(3) 



When we evolve to high Q 2 we probe the proton 
structure ever more finely, to transverse sizes ~ W Q 2 , 
see Fig. 2. 

The parton distributions are essential to calculate 
the cross sections, a, for "hard" hadronic processes. 
First the QCD subprocess, a, are calculated in the 
strongly-ordered kj, = approximation and then the 
factorization theorem gives a = xg(x,/j, 2 ) ® ^{l^ 2 ) + ■■• 
where the scale fj, 2 ~ hard scattering pj,. 

(b) BFKL equation (small x) 



(a) GLAP evolution (large Q 2 ) 

Starting from a known structure of the proton at 
Q 2 = Qo we ma y evolve up to large logQ 2 using the 

* To be published in the Proc. of ICHEP Conf., Glasgow, July 
1994. 



Figure 1. A sample of the latest HERA data for F 2 [1,2]. The 
GLAP-based MRS(A) analysis [3] and the BFKL-based AKMS 
"prediction" [4] give almost indistinguishable descriptions of the 
HERA data. Also shown are the GRV parton [5] and "hot-spot" 
shadowing [6] predictions. 
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Figure 2. The gluonic content of the proton as "seen" in 
different deep-inelastic (x, Q 2 ) regions. W is the ratio of the 
quadratic to the linear term on the right hand side of (10). 



On the other hand, when we evolve up to large 
1/x (i.e. small x) we encounter (a s log(l/x))" terms 
which have to be resummed. Indeed the dramatic rise 
observed in F2 with decreasing x may be associated 
with the growth of the gluon density which arises from 
the resummation of these terms; a growth which, via 
g — > qq, is transmitted to the sea quarks probed 
by the photon. To leading order the summation is 
accomplished by the BFKL (or Lipatov) equation ||, 
which may be written in the differential form 
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There is now no strong-ordering in k n T of the emitted 
gluons and we have to work in terms of the unintegrated 
gluon distribution f(x,k T ) in which the "last" k T 
integration along the gluon ladder (of Fig. 3a) is 
unfolded, that is 



n* 2 dfc 2 

xg(x,fi 2 )^J -^ff(x,k T ). 



(5) 



At small x the gluon distribution f(x,k T ,ji 2 ) becomes 
independent of the scale fx 2 . 

From (4) we see that the small x behaviour of 
/ is controlled by the largest eigenvalue Xl of the 
eigenf unction equation K ® f n — A„/„, since as x —* 

/~exp(A L log(l/x))~x- Ai . 

Indeed for fixed a s there is an analytic solution for the 
leading small x behaviour 



/ - x-* L (k T )*exp -c 



\og 2 (k 2 T /k 2 ) 



where 



log(l/x) 
A L = (3a s /7r)41og2. 



(6) 
(7) 



Figure 3. (a) Gluon ladder, (b) diagrammatic representation of 
the BFKL contribution to F2, see (13). 



This singular x L Lipatov behaviour is in contrast to 
the naive Regge-type expectations that 
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A-a P (0) 



-0.08 



(8) 



where ap(0) is the intercept of the Pomeron. 

A second feature of the solution (6) of the BFKL 
equation is the diffusion in kx with decreasing x, as 
manifested by the Gaussian form in logfc^ with a width 
which grows as (log(l/x))^ as x decreases. The physical 
origin of the diffusion is clear. Since there is no strong- 
ordering in kr, there is a "random walk" in kr as we 
proceed along the gluon chain and hence evolution to 
smaller x is accompanied by diffusion in kx- We foresee 
that the diffusion will be a problem in the applicability 
of the BFKL equation since, with decreasing x, it 
leads to an increasingly important contribution from the 
infrared and ultraviolet regions of k\ where the equation 
is not expected to be valid. 

For running a s the singular behaviour and diffusion 
in kx are confirmed. In addition it is found that 0] 



/ ~ C(k T )x 



-A 



(9) 



where the value A w 0.5 is much less sensitive to the 
treatment of the infrared region in (4) than is the 
normalization C. 



(c) Shadowing region 

The x~ x growth of the gluon cannot go on 
indefinitely with decreasing x. It would violate unitarity. 
The growth must eventually be suppressed by gluon 
recombination, which is represented by an additional 
quadratic term so that (4) has the symbolic form 



xdf/dx = K(g>f-V(g)f 2 



(10) 



The additional term contains a factor a 2 /k^R 2 since the 
gluon-gluon interaction behaves ~ a 2 /k T , whereas 1/R 2 
arises because the smaller the transverse area (ttR 2 ), in 
which the gluons are concentrated within the proton, the 
stronger the effect of recombination. The precise form 
of this equation, originally proposed by GLR ||, is still 
a matter of debate (l(J. The region where shadowing 
should be calculable perturbatively is just below the 
dashed line in Fig. 2. 

2. Small x behaviour of Fi 

The small x behaviour of xg (and xq) arising from 
GLAP evolution depends on the form of the starting 
distributions. For singular starting distributions, xg ~ 
x~ x with A > 0, the small x behaviour is stable to 
evolution in Q 2 . The larger the value of A the sooner the 
stability sets in with decreasing x. MRS (A) partons S, 



with xg 



-0.3 



, are an example of this behaviour. On 



the other hand, for non-singular starting distributions, 
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xg ~ : 
logarithm (DLL) form 



with A < 0, we find the double leading 



a; ff ~exp(2[£(Qg,Q 2 )log(i/x)]5). 



(11) 



That is xg grows as x — > faster than any power of 
\og(l/x) but slower than any power of x. The larger the 
"evolution length" , 



Q2 df 3a s (q 2 ) 
Ql I 2 * 



(12) 



the faster the growth. An example is the "dynamical" 
GRV partons || which evolve from valence-like forms at 
a low scale Q\ — 0.3 GeV 2 , and for which (11) mimics 



a behaviour xg 



-0.4 



in the HERA regime. 



Given the solution f(x,k T ) of the BFKL equation 
we can use the fcT-factorization theorem to predict F 2 , 
see Fig. 3b: 



F 2 = f®F h 
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(13) 



where A ~ 0.5, and F 2 6 ~ 0.4 is determined from the 
large x behaviour of F%. Once the overall normalisation 
of the BFKL term is adjusted by a suitable choice of the 
infrared parameters in (4), then an excellent description 
of all the F 2 (x,Q 2 ) HERA data is obtained. Indeed 
the BFKL-based "prediction" j| gives an equally good, 
and almost indistinguishable, description as the GLAP- 
based order fit ||, see Fig. 1. With GLAP, the steepness 
is either incorporated (as a factor x~ x ) in the starting 
distributions or generated by evolution from a low scale 
Qq. The steepness can be adjusted to agree with the 
data by varying A or Q 2 . On the other hand the leading 
log(l/x) BFKL prediction for the shape F 2 — i 7 ^ 8 ~ x~ x , 
with A ~ 0.5, is prescribed. It remains to see how well 
it survives a full treatment of sub- leading effects. 

Conventional shadowing with gluons spread uni- 
formly across the proton (R — 5 GeV -1 ) leads to only a 
small suppression in F 2 in the HERA regime. If the glu- 
ons were concentrated in "hot spots" of area ttR 2 with, 
say, R = 2 GeV -1 the effect would be much stronger, 
see Fig. 1. But could shadowing be identified since we 
do not know the partons at small xl Simulated F 2 data 
(of accuracy and x range which may eventually be ac- 
cessible at HERA) have been used jllj to see how well 
R could be determined. The conclusion is that the in- 
terplay between the linear and non-linear terms in (10) 
leads to a considerable ambiguity between the size of 
the parton distributions and the amount of shadowing. 

Is GLAP evolution adequate in the HERA regime? 
For sufficiently small x the (a s logl/a;) n terms must be 
resummed with the full Q 2 dependence (and not just 
the leading and next-to- leading logQ 2 terms). Ellis et 
al. |Q have made a theoretical study of the applicability 
of GLAP evolution and find that it is adequate in 



Figure 4. Processes that may be used to identify BFKL 
dynamics. 



the HERA small x regime, provided that the evolution 
occurs from a sufficiently singular starting distribution, 
xg ~ x~ x with A > 0.35. 



3. Identification of BFKL behaviour 

The inclusive nature of F 2 , and the necessity to provide 
"non-perturbative" input distributions of parton densi- 
ties for its description, prevents its observed small x be- 
haviour being a sensitive discriminator between BFKL 
and conventional dynamics. For this purpose it is nec- 
essary to look into the properties of the final state. 

The two characteristic features of BFKL dynamics 
are the absence of strong-ordering of the gluon k^s 
along the chain (the diffusion in kx) and the consequent 
(x/x')~* or exp(AAy) growth of the cross section, where 
x and x' are the longitudinal momentum fractions of 
the gluons at the ends of the chain, which spans the 
rapidity interval Ay = log(a:'/a;). Recall A ~ 0.5. Some 
processes which exploit these characteristic features are 
shown in Fig. 4. 

The idea [jl3| in Fig. 4a is to detect deep-inelastic 
(x,Q 2 ) events which contain a measured jet (xj, k T j) 
in the kinematic regime where (i) the jet longitudinal 
momentum, Xj, is as large as is experimentally feasible 
(xj ~ 0.1), (ii) z = x/xj is small, and (iii) k T j ~ Q 2 is 
sufficient to suppress diffusion into the infrared region. 
The beauty of this measurement is that attention is 
focussed directly on the BFKL z~ x behaviour at small z. 
The difficulty is to cleanly separate the forward going jet 
from the proton remnants. The preliminary results from 
the HI collaboration jl4| are encouraging and favour the 
BFKL over the conventional description. 

Inspection of Fig. 4b suggests, that due to the 
relaxation of strong-ordering of the gluon fc^'s at small 
x, more transverse energy Et should be emitted in 
the central region (between the current jet and the 
proton remnants) than would result from conventional 
evolution. Indeed Monte Carlo predictions based on 
QCD cascade models]]] fall well below the observed 
central plateau of height Et ~ 2.1 GeV per unit of 
rapidity |ll|. A BFKL-based calculation Q, at the 
parton level, yields about 1.7 GeV per unit of rapidity, 
but much less if conventional dynamics is used. No 
hadronization effects have been allowed for. 

f A good description has been obtained by a Monte Carlo based 
on the colour dipole model. This model contains the essence of 
BFKL dynamics provided the full integration over the emitted 
gluon kx is performed. 
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Recently there has been renewed interest in the 
original proposal of Mueller and Navelet Jl7| that the 
cross section for the production of a pair of minijets 
should, according to BFKL dynamics, rise as exp(AAy) 
as the rapidity interval Ay becomes large. The studies 
jiJI show that the effect is masked by the fall-off of the 
parton densities at large x, but that instead the rate of 
weakening of the azimuthal (back-to-back) correlation 
between the jets, could possibly be an indicator of BFKL 
effects. 

BFKL dynamics may be also identified via the 
weakening of the azimuthal correlation between a pair 
of jets produced in deep-inelastic scattering at HERA, 
see Fig. 4d. At sufficiently large values of A<fi = 
<j) — 7T, BFKL dynamics dominates over the fixed- 
order QCD contribution from 3+1 jet production, 
leading to a distinctive tail in the azimuthal distribution 
which directly probes the kr dependence of the gluon 
distribution Q. 

4. Conclusions 

In the HERA small x regime GLAP evolution (from 
appropriately parametrized starting distributions) is 
able to mimic BFKL dynamics as far as the description 
of F2(x,Q 2 ) is concerned. Moreover it will be difficult 
to isolate shadowing contributions even with improved 
measurements of F%. Measurements which are less 
inclusive than F2, offer more chance to identify the 
characteristic BFKL x~ x behaviour and diffusion in kx- 
However opening up the final state brings the problems 
of hadronization and jet identification, and loses some of 
the small x "reach" (e.g. x — > x/xj where xj ~ 0.1). On 
the theoretical side, the sub-leading corrections to the 
BFKL leading log(l/:r) formalism are urgently needed 
for future quantitative studies of small x phenomena. 
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